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Ph! Abstract 

<D . 

' The extension of strongly anisotropic or dynamical scaling to local scale 

invariance is investigated. For the special case of an anisotropy or dynamical 

* 1 

exponent 9 = z = 2, the group of local scale transformation considered is 
^ the Schrodinger group, which can be obtained as the non-relativistic limit of 

the conformal group. The requirement of Schrodinger invariance determines 
the two-point function in the bulk and reduces the three-point function to a 
scaling form of a single variable. Scaling forms are also derived for the two-point 
function close to a free surface which can be either space-like or time-like. These 
results are reproduced in several exactly solvable statistical systems, namely 
the kinetic Ising model with Glauber dynamics, lattice diffusion, Lifshitz points 
in the spherical model and critical dynamics of the spherical model with a non- 
conserved order parameter. 

For generic values of 0, evidence from higher order Lifshitz points in the 
spherical model and from directed percolation suggests a simple scaling form 
of the two-point function. 



present adress 



1 Introduction 



Scale invariance is a central notion in present theories of critical behaviour. In the 
context of two-dimensional, static and isotropic critical behaviour, these ideas have 
become spectactularly succesful in the context of conformal invariance [p]] . The main 
physical idea behind this is the extension of the covariance of correlation functions 
under length rescaling by a constant factor A to general, space-dependent rescalings 
A(r ). A coordinate transformation A(r ) is conformal if the angles are kept unchanged. 
In two dimensions, the conformal group is infinite-dimensional. This has lead, for 
example, to the exact calculation of critical exponents and correlation functions and 
yields a handle for classifying two-dimensional universality classes (for reviews, see 



Much less is known about non-isotropic scaling. Consider a (connected) correla- 
tion function C(f;t) depending on "space" coordinates r and a "time" coordinate t 
which satisfies the scaling relation 



where x is a scaling dimension and 9 is referred to as anisotropy exponent. Systems 



terns. In fact, dynamical scaling of this kind apears quite commonly in time-delayed 
averages close to an equilibrium phase transition, where the anisotropy exponent is 
referred to as dynamical exponent z = 8, see fg, [|], or in domain growth problems of 
systems quenched below the equilibrium critical point, see [||, || and for a recent 
review. Alternatively, strongly anisotropic scaling may arise in statics, examples are 
provided by directed percolation || or by magnetic systems at a Lifshitz point 0, 
where the anisotropy exponent 9 = v\\/v±_ is related to the critical exponents u»± 
of the correlation lengths £im parallel and perpendicular to the preferred direction. 
Eq. (|1 . 1|) can be recast in the form 





C{Xf;\ e t) = \~ 2x C{f-t) 




which satisfy eq. ( |1 . 1[ ) with 6^1 are by definition strongly anisotropic critical sys 
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which defines the scaling function and 

u = j (1.3) 

is the scaling variable and henceforth we shall be always taking the scaling limit 
r — > oo, t — > oo where u is kept fixed. 

We ask the following question: what can be said about the scaling function $(w)? 
Is it sensible to look beyond global scaling with A constant to a space-time dependent 
rescaling factor X(r,t)7 

This question had been adressed by Cardy ||10|| . Assuming dynamical scaling for 
the dynamic response function, he takes as extended set of local scale transformations 
A(r) the space-dependent scaling r — > A(r)r, t — > X(f) e t where the A(r) are two- 
dimensional conformal transformations. This means that only systems at a static 
critical point are considered. The assumed covariance of the response function is used 
to map the problem from the two-dimensional plane (in the space coordinates) to the 
strip geometry, with a non-uniform rate. Next, since close to a static critical point the 
static correlation length of the system in the strip is of the same order of magnitude 
as its width, it is claimed that . . on much larger distances it is permissible to use 
mean field theory to calculate the dynamic correlation function in the strip" [TO]. For 
a system with a non-conserved order parameter the response function then turns out 
to be [|Kj 

G(r, t) ~ t - 2x/e - 1 exp (-jj (1.4) 

where some non-universal constants have been suppressed. The case of a conserved 
order parameter was also treated. While this result is appealingly simple, the as- 
sumptions made in deriving it may appear to be quite strong,^ in particular the use 
of mean-field (van Hove) theory. Also, one might wish to reconsider the assumption 
that A(r ) were time-independent. In fact, we shall study the scaling of the two-point 
function in (1 + 1)D directed percolation and higher order Lifshitz points in the 
2 The restriction to 2D space is not really required and could be removed 
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spherical model, both with 8^2, and find that the form of the scaling function of 
these models does not agree with eq. (|1.4j) . 

Here we propose another group of space-time dependent local scalings A(r, t). For 
definiteness, we shall only consider the case = 2, but we do not have to make the 
restriction to d = 2 space dimensions. Although the van Hove theory for a non- 



conserved order parameter has also 6 = 2 (see [12fl), we do not make any of the 
approximations involved in that theory. Rather, it is our aim to find the scaling func- 
tions merely from their transformation properties under local scale transformations. 
By taking 9 = 2, we mean to perform the simplest case study of local, albeit not con- 
formal, scaling transformations. The group of the local scaling transformations is the 
Schrodinger group, which shall be defined in the next section. The approach chosen 
has the advantage of being close in spirit to the earliest investigations of conformal 



invariance in critical phenomena [ 13| . On a more formal level, the comparison of 



conformal with Schrodinger invariance provides some insight into the characteristic 
properties of both. 

Our results are as follows: 

1 . If the domain of both time and space coordinates is infinite in extent and the 
scaling fields transform covariantly under the Schrodinger group, the two-point 
function is completely determined, while the three-point function is reduced to 
a scaling form of one variable, see eqs. ( |3.12| , |3.28D . For example, this applies to 
the calculation of time-delayed correlations of systems at equilibrium and at a 
static critical point, or else to lattice diffusion problems. 

2. If the space geometry is semi-infinite, a scaling form for the two-point function 
will be derived, see eq. (|3.35| ). This may be relevant to critical dynamics close 



to a surface, see for an example. 



3. For a system being in a predefined initial state, it can be shown that critical 
relaxation towards equilibrium displays dynamical scaling already at interme- 
diate times 0, much later than microscopic times but also well before the late 



time regime usually considered. We derive a scaling form for the two-point 
function, see eq. ( |3.44| ). 

4. These results can be reproduced from a variety of exactly solvable models. 

5. The Schrodinger group can be extended to an infinite-dimensional group, whose 
Lie algebra contains a Virasoro subalgebra. It can be shown that for systems 
with local interactions Schrodinger invariance follows from the requirements of 
translation invariance in both space and time, rotation invariance in space, scale 
invariance and Galilei invariance. If no anomalies occur, this even holds for the 
whole infinite-dimensional group. 

The work described in this paper uses background from both conformal field 
theory and time-dependent statistical mechanics. To make the paper accessible 
to readers with knowledge in one, but not both of these fields, we repeat in sec- 
tion 2 the definition of the Schrodinger group and recall a few well-known facts about 
Galilei-invariant theories and dynamical scaling. Section 3 describes the derivation 
of the two- and three-point functions for either infinite or semi-infinite geometries. 
The Schrodinger Ward identity is considered as well, and the non-existence of non- 
conventional central extensions of the Schrodinger Lie algebra is shown (Appendix B). 
The discussions of this section follow closely the known derivation of correlation func- 
tions from conformal invariance. In section 4, we test and confirm the predictions 
from Schrodinger invariance by calculating two- and three-point functions in several 
exactly solvable and strongly anisotropic critical models. In section 5, we examine 
the scaling of the two-point function for some systems with an anisotropy exponent 
6 2. Exact and numerical results indicate a disagreement with that result eq. ( |1.4|) 



obtained from 2D conformal invariance and suggest an alternative simple scaling 
form. Section 6 gives our conclusions. 
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2 Basic concepts and terminology 



We begin by recalling some well-known facts about the Schrodinger group, Galilean 
invariance in field theories and dynamical scaling. 

2.1 The Schrodinger group 



The Schrodinger group is defined |15| [L6[] by the following set of transformations 

izf+vt + d 



, at + P 
t — > t — — : ; ad — p7 = 1 



(2.1) 



jt + 5 -ft + 5 

where a, j3, 7, 5, v, a are real parameters and TZ is a rotation matrix in d space dimen- 
sions. It is apparent that the Schrodinger group can be obtained from the Galilei 
group by extending the time translations to the full Mobius group Sl(2,R) of frac- 
tional real linear transformations in time as given in eq. (|2.1| ). A faithful matrix 
representation is given by 



v a 



( U 

a /3 

\ 7 5 



(2.2) 



Niederer [yj| showed that this group is the maximal kinematical group which 
transforms solutions of the free Schrodinger equation 



(i— + ——) ib - 
y dt 2m dr 2 J 

into other solutions of (|2.3| ), viz. (f,t) > g(r,t), tb T g ib 

(T g ib)(r,t) = f 9 (g~\r,t)) ib (g-\?,tj) 
where the companion function f g is [|l^ 



(2.3) 



(2.4) 



/ 9 (f,t) = ( 7 t + 5)- d/2 exp 



im 7r 2 + 21Zr ■ (7a — 5v ) + 7a 2 — 5tv 2 + 2^av 



-ft + 5 



(2.5) 



Independently, it was shown by Hagen ]T6|] that non-relativistic free field theory is 
Schrodinger invariant, treating both scalar and spin-^ fields. It was also shown that 



the operators which appear in the conservations laws associated with the space-time 
symmetries can be reformulated to allow the statement of Schrodinger invariance 
entirely in terms of those operator densities The Schrodinger group in d space 
dimensions can be obtained by a group contraction (where the speed of light c — > oo) 
from the conformal group in d + 1 dimensions ]T7| provided the mass is conviently 
rescaled as well. Its projective representations as required by (|2.4j) have been studied 



in detail [IS 



There are many more equations, whose kinematical group is isomorphic to the 
Schrodinger group. It can be shown that the most general potential which can be 
added in ( |2.3| ) such that the kinematical group is still isomorphic to ( |2.1| ) is of the 
form, up to orthogonal transformations V(r) = J2i r f + Hi v i r i + v ^ 101 ■ ^ur- 
ther examples are provided by a non-linear Schrodinger equation or, but with a more 
general transformation law than fl2.4|) , by the Navier-Stokes equation with homoge- 
nous pressure or by Burger's equation [pOfl . Higher order symmetry operators of 



eq. ( |2.3| ) are examined in |2T|. The Schrodinger group also appears as a dynamical 
symmetry group for the Dirac monopole or magnetic vortices or in the non- 



relativistic N— body problem with inverse-square interactions [23]. Similarly, one 
may treat the diffusion equation by writing m _1 = 2iD, where D is the diffusion con- 
stant. We shall do so for most of this paper. In any case, we shall only consider here 
the realization of the Schrodinger group provided by the free Schrodinger equation 

!■ 

For simplicity, we restrict attention here to fields which are scalar under space 
rotations (so that it is sufficient to take 1Z — 1) and shall mostly also take just 
one space dimension d — 1. This is not a serious restriction and generalizations are 
straightforward. 

The set Sa n = {X_i, X , X\, YL1/2, Y1/2, M } spans the Lie algebra of the Schro- 
dinger group eq. ( |2.1| ). The generators read (we take d — 1) 

Xn = -t n+1 d t - ^±t n rd r - ^ + l) Mt n ~ l r 2 

2 4 
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Y m = -t m+1/2 d r - (m + A4i m " 1/2 r 

M n = -t n M (2.6) 

where the terms ~ M. come from the companion function. When M. = im is purely 
imaginary, this corresponds to the Schrodinger equation where m is the mass, while 
for Ai real, this is the form corresponding to the diffusion equation. The commutation 
relations are 



[X n ,X m ] = (n-m)X, 



n+m 



n 



\X n , Y m \ i ^ ^ TTij Y n j rm 

[X n , M m ] = -mM n+m 

[Y n , Y m ] = (n- m)M n+m 

[Y n ,M m ] = [M n ,M m ] = (2.7) 

(In more than one space dimensions, there are several sets of generators Y$ , i = 
1, . . . , d, but only one set of X n , where rd r is replaced by rd? etc.) The commutation 
relations ( |2.7| ) remain valid when the infinite set of generators S = {X n , Y m , M n } 
where n is an integer and m is a half- integer, is considered The Lie algebra can 
be decomposed S = Sx ® s <Sy where Sx = {X n } and Sy = {Y m , M n }. As we shall 
see later, these two subalgebras arise in quite distinct physical situations. 

2.2 Galilei-invariant field theory 

The Galilei Lie algebra (here for the case d = 1 only) is generated from the set 
Q = {X_x, y_i/ 2 , Y1/2, } C Sa n . Thus any Schrodinger invariant theory will have to 
satisfy the constraints following from Galilean invariance as well. These conditions are 



well known for a long time |25] and we briefly recall the properites relevant for us. In 
fact, it is possible to construct a consistent field theory which is Galilei-covariant from 
the following postulates. The simplest example of this is second-quantized ordinary 
non-relativistic quantum mechanics. 
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States are rays in a Hilbert space and the dynamical variables are operators in the 
Hilbert space, with the usual rules for the calculation of probabilities. The Galilei 
group acts by a unitary projective representation U(g) in the Hilbert space. If <p(r, t) 
is a field of the theory, it is required to transform locally 



U(g)- 1 <l>(r,t)U(g) = exp 



vi i v , n ~ \ 

y(tit + 2vr) 



{r + vt + a,t + (3) (2.8) 



States are characterized by the Casimir operators, whose eigenvalues are mass and 
spin (if d > 1). Since we deal with projective representations, we obtain a unitary 
representation of a central extension of Q by a one-dimensional Lie algebra generated 



by M , see eq. (|2.7| ). This implies, since the extension is non-trivial because Q is 
non-semisimple, that a physically trivial transformation may result in a modification 
of the phase of the state vector which depends on the mass of the system. Galilei 



covariance thus requires the Bargmann superselection rule of the mass |26| which 
states that for an interaction of particles of the form 

A + B + > A' + B' + ■■■ (2.9) 

one must have 

m A + m B -\ = m A ' + mB'-\ (2.10) 

This implies that no Galilean field can be hermitian unless it it massless. We see 
that the mass plays quite a distinct role in non-relativistic theories as compared to 
relativistic ones. We emphasize that the mass no longer describes the deviation from 
critical behaviour in our context, as it does in relativistic theories. Masses should, 
in the light of ( 2.10Q , rather be regarded as some kind of analogon of a charge [25 |. 



Non-vanishing correlations are of the type 

{cf> a {r,t)<j>t{r',t')) ~ 5„ b X a , b (r- f',t- t') (2.11) 

and similarly for higher order correlations. This also holds by analytic continuation 
for Euclidean theories. We shall rederive the superselection rule eq. ( |2.10|) in several 
cases below. 
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A further remark is in order here. In general, in the context of a statistical 
system, the mass will contain non-universal factors which merely serve to define the 
time scale. Here we are interested in the ratios of masses of different scaling fields, 
which are universal. 

While the Bargmann superselection rule provides a restriction not present in rel- 
ativistic theories, Galilean field theory is considerably less restricted in many other 
aspects. For example, the consequences of locality in Galilean field theory are no 
longer sufficient to prove neither the CPT nor the spin-statistics theorem, see for 
a full discussion. 

2.3 Dynamical scaling 

Finally, we recall some facts about response functions and dynamical scaling, follow- 



ing ][|, 27]. Consider a scaling field (f)(r,t). In systems described by a Hamiltonian 
one can define its conjugate field h(f,t). Then the linear response function x(k,u) is 
defined in momentum-frequency space by 

< <f)(k, to) >= x(k, u)h(k, uj) (2.12) 

where the field h is taken to be infinitesimal, the average is determined from the time- 
dependent probability distribution in the presence of h and the system is assumed to 
start from equilibrium at t — > — oo. The Fourier transforms are 

X{k,uj) = J o dt J dfe i{ujt - k ^G(r,t) (2.14) 
We are interested in time-delayed correlation functions 

CV(f, t) =< 0(f, t)0(O, 0) > h=0 - < 0(f, t) > h=0 < 0(0, 0) > h=0 (2.15) 



and define its Fourier transform C^(k,u) according to eq. ( |2.13|) and the equal-time 



correlation function C^k) as the spatial Fourier transform of (2.15|) with t = or 



Cm= — C^(A;, a;) (2.16) 

J-oo Z7T 



From causility, it can be shown (see [27]) that the response function x(k,oj) is an 



analytic function of the complex frequency u in the upper half-plane and its real 
and imaginary parts satisfy the Kramers- Kronig dispersion relations. For classical 
systems with a Hamiltonian, the fluctuation-dissipation theorem states [§], ^7]] 

C,(£ )W ) = ^_3x(*,w) (2.17) 

where 3 denotes the imaginary part. The hypothesis of dynamical scaling now asserts 
that at a static critical point |3|, [| 

X (k, u) = A k 2x - d §(Bujk- d ) (2.18) 

in the scaling limit uo — > 0, k — > with uk~ e fixed, where 9 = z is the dynamical 
(anisotropy) exponent, x a scaling dimension, $ is a universal scaling function and 
A, B are non-universal constants. 

3 Multipoint correlations from Schrodinger invari- 
ance 

We now turn to derive the consequences of Schrodinger invariance for the correlations. 
In general, we expect a scaling field 4>(r,t) to be characterized by its mass Ai, its 
scaling dimension x and its spin s (which we take to be zero throughout, but see 
16H for the case of spin 1/2). The discussion will be exclusively for d — 1, but the 
extension to arbitrary d is immediate. The transformation of <f)(r, t) will contain 
terms describing the space-time coordinate change given by A(r, t), the scaling as 
described by the Jacobian of A(r, t) and the change in the phase which is a peculiar 
feature of non-relativistic systems. Under infinitesimal coordinate changes, we have 
the transformations 

[X n ,<f>(r,t)] = [t n+1 d t + ^t n rd r + ^W -V + (n + <P(r,t) 

\y m ,4(r,t)] = (^t m+1 /%+(m+^Mt m - 1 / 2 r^ ( j ) (r 7 t) (3.1) 
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Taking over the conformal terminology of 0, we call a field primary, if it satisfies ( |3.1|) 
for all n integer and all m half-integer. A field is called quasi-primary, if it satisfies 
(|3.1D for the finite-dimensional subalgebra Sa n only. Consider multipoint correlators 
{^ai^a, ~ta)4>b( r bi h) ■ ■ ■ 0y( r ^) ty)4 ) * z ( r z, t z )) of quasiprimary fields and we derive the re- 
strictions following from the hypothesis of their covariant transformation under 5g n . 
We shall use the short-hand 

(3 ' 2) 

We shall not consider explicitly the action of the generator Mo, because invariance 
with respect to it follows from the Bargmann superselection rule. 

3.1 Two-point function in the bulk 

We consider the two-point function 

F = F(r a , r b ; t a , t b ) = (<p a (r a , t a )(p* b (r b , t b )} (3.3) 

of quasiprimary fields 0^ in the infinite geometry in both time and space. Invariance 
under translations in time and space implies F = F(r, t), where r = r a —r b , r = t a —t b . 
Invariance under scale transformations generated by X requires 

(tad a + \r a D a + ^ + t b d b + X -r b D b + ^) F(r, r) = (3.4) 

which is rewitten as, with x = \{x a + x b ) 

(rd T + -rd r + x) F(r, r) = (3.5) 

We write the solution in the form 



^2 



F(r,T) = T- x G\^—J (3.6) 

which is nothing but the scaling form eq. ( |1.2| ). New information comes from requiring 
Galilei invariance (V1/2) 

(t a D a + M a r a + t b D b - M b r b ) F(r, r) 
= (rd r + M a r a -M b r b )F(r,r) = (3.7) 
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and we obtain two conditions 

f M a -M b = 



(rd r + M a r) F(r, r) = ( 3- * 



We recognize in the first of these the Bargmann superselection rule eq. (|2.10| ). Com- 
bining with scale invariance (|3.6| ) we find 



G(u) = G exp (-^u) (3-9) 

We remark that the form of F(r, r) in ( |3.6|) is not an arbitrary ansatz. This can be 
seen by first solving the conditon of Galilei invariance before using scale invariance. 
Finally, invariance under the special Schrodinger transformation X\ gives 

(f a d a + t a r a D a + M a r 2 a + x a t a + t 2 b d b + t b r b D b - M b r 2 b + x b t b ) F(r, r) = (3.10) 

which is seen as before to lead to the conditions 

M a -M b = (3.11) 
G' + \M a G = 

where the prime denotes the derivative. The last two of these had already been 
obtained before. The final result is, where $o is a normalization constant 

M a (r a - r b f 



F = S XaiXb 6 MaMb ^ (t a - t b ) Xa exp 



(3.12) 



2 t a — t b 

which should be understood in the scaling limit. This had been announced before 



24}| for the special case of equal masses. Since for a non-conserved order parameter 
van Hove theory leads to a dynamic exponent z = 2, it is not surprising that we 
recover in this case the form eq. ( |1.4j ) as found for d = 2 by conformal invariance [IU|. 



3.2 Three-point function in the bulk 

Consider the three-point function 

F = F(r a , r b , r c ; t a , t b , t c ) = (<fr a (r a , t a )(p b (r b , t b )(p* c (r c , t c )) (3.13) 
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Translation invariance in both space and time let F only depend on distances, F = 
F(r, s;r,a) where 

r = r a -r c , s = r b - r c , r = t a -t c , a = t b - t c (3.14) 

Scale invariance requires that 

Yl (tidi + 2 TiDi + 2 Xi ) F ^ T ' S] T ' ^ 
= (rd T + ad a + ^rd r + ^sd s + - (x a + x b + x c )\ F = (3.15) 

Making the ansatz 

F(r,s;r,a) = T-' J1 a- p2 (r - a)- p[i G(r,s;r,a) (3.16) 

as motivated by the corresponding result for the three-point function as obtained 
from conformal invariance we find 

j Pl+ P2 + Pz = \{Xa + Xb + X c ) 

\ (rd T + od a + \rd r + lsd s ) G(r, s; r, a) = 1 ' 

Note that the relation found between the exponents pi and the scaling dimensions Xi 
must be satisfied for any scale-invariant system. From Galilei invariance we get, with 
e a = e b = —e c = 1 because of eq. 



5^ (*iA + aMiTi) F(r, s; r, a) (3.18) 

i=a 

= (rd r + ad s + M a r + M b s + (M a + M b - M c ) r c ) F(r, s; r, a) = 

which leads to the conditions 

f M a + M b -M c = , . 

\ (r<9 r + a9 s + .M a r + >f 6 s) G(r, s; r, a) = l ^ Jj 

and we recognize again the Bargmann superselection rule for the masses. The equa- 
tion for G can be further simplified by setting 

G(r, s; r, a) = exp " H ^ s ' r ' a ) ( 3 ' 2 °) 
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Since the first factor is scale invariant, H satisfies the same eq. ([3.171) as does G, but 
the second eq. ( |3.19| ) becomes 



(rd r + ad s ) H(r, s; r, a) = 



(3.21) 



Finally, invariance under the special Schrddinger transformation requires 

(tidi + UnDi + -eMir* + t iXt ) F(r, s; r, <r) = (3.22) 

i=a - / 

First, we use the form ( |3.16 ). Then we get (T + T>)G = 0, where T> is a differential 
operator and 



T = t a {-pi - p 3 + x a ) + t b (-p 2 - p 3 + x b ) + t c {-p x - p 2 + x c 



(3.23) 



The requirement that T vanishes leads together with scale invariance (|3.17|) to 

P2 = ^(x b + x c - x a ) , pi = ^(x a + x c - x b ) , p 3 = ^(x a + x b - x c ) (3.24) 

Secondly, we use eq. (|3.20|) and get (T 7 + V')H = where 

1 



(3.25) 



This vanishes due to the Bargmann superselection rule. The differential operator V 
finally takes the form where we used that H is scale as well as Galilei invariant 



(r 2 d T + a 2 d a + rrd r + <Jsd s ^ H(r, s; r, a) = 



(3.26) 



It remains to solve the resulting system of first order linear partial differential equa- 
tions. This is done in Appendix A, with the result 



H(r,s;r,a) = ^ 



(ra — st) 



(3.27) 



[a — r )ar J 

where \1/ = ^ a b,c is an arbitrary function. The final result is, with the appropriate 
scaling limits understood 



F = 5 



Ma+M b ,M 



(t a ~ t c )-^ Xa+Xc ~ Xb \t b - t c )-^ Xb+Xc - Xa) (t a - t b )-^ Xa 



+x b -x c ) 



x exp 



M a {r a - r c f M b (r b - r c f 



2 th — t r 



'[{r a ~ r c )(t b - t c ) - (r b - r c )(t a - t c )] : 



\pa — tb)(ta — t c )(t b — t c ) 



(3.28) 
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In particular, it follows that any three-point function < 000* > of a massive field 
with itself vanishes. 

The results eqs. ( |3.12| , |3^28D deserve some comments. It is instructive to compare 



them with those for the two-point and the three-point functions obtained from the 
requirement of conformal invariance by Polyakov see also |2^|. Conformal 
invariance completely specifies the form of the two- and three-point functions in any 
number of space dimensions. Also, two-point correlations of quasiprimary fields must 
vanish if the scaling dimensions are different. We reproduce this result in eq. (|3.12| ), 



but add the stronger requirement of the Bargmann mass selection rule. The exponen- 
tial behaviour of the two-point scaling function is a consequence of Galilei invariance, 
while $o is merely a normalization. Our results do depend on the explicit realization 
of the Galilei transformation as given by the generator Yxi%- If we had considered 
the Schrodinger group with a potential present, the realization of the generators is 



different [19) and we would have found a different form of the correlations. Turning 



to the three-point function eq. fl3.28| ), we observe that the purely time-dependent 



factors reproduce the familar form of the conformal three-point function which is 
completely symmetric in the times t a ,ti,,t c . We then note that this symmetry is not 
met by the exponential factors, determined from Galilei invariance and we also note 
the appearance of the Bargmann mass selection rule, not present in conformal invari- 
ance. Generalization of our results to higher space dimension d merely requires to 
check for rotation invariance, which is obviously satisfied. 

3.3 Two-point function in semi-infinite space 

Having studied correlations of quasiprimary fields in the infinite geometry, we now 
consider the effect of surfaces. Consider a free surface at r = 0. It is kept invariant 
under the transformations generated by the subalgebra Sx, but space translations 
and Galilei transformations will no longer leave the system invariant. Nevertheless, it 
is known that conformal invariance can be used in analogous situations to constrain 



the two-point correlation function p9fl . For quasiprimary fields, we require covariance 
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only under the subalgebra Sx H <Sfi n = X , Xi}. 

Consider the two-point function of quasiprimary fields 



F = F(r a , r b \ t a , t b ) = (0 a (r o , t a )(f)l(r b , t b )) 



(3.29) 



and we require space points to be in the right half-plane, i.e. r a , r b > 0. Time 
translation invariance gives F = F(r a , r b ; r), with r = t a — t b . From scale invariance 
we obtain 



£ (tA + ir<A + ~ Xi ) F 
r9 T + X -r a D a + ^7,/^, + ./■ j F = 



(3.30) 



where x = |(x a + x b ). On the other hand, from the invariance under the special 
Schrodinger transformation we have, with e a = —e b = 1 

J2 {t 2 A + UnDi + l -e l M % r] + Ux^j F 

T 2 d T + rr a D a + t b (2rd T + r a D a + r b D b ) 
+ \M a r 2 a - ^M b r 2 b + t a x a + t b x b ^j F 
T 2 8 T + rr a D a + l -M a r 2 a - l -M b r 2 + rx^ F = 



(3.31) 



where in the last equation the scale invariance of F was used. Now, we make the 
ansatz 



F(r a ,r b ;r) = r x G(u,v) , u 



" v = ^ 



T 



which solves for scale invariance, while eq. (|3.31|) gives 



X x a x b 
(ud u - vd v + \M a u - \M b v^j G(u, v ) = 



(3.32) 



(3.33) 



The general solution of this is found using the method of characteristics [3T)| 



G(u, v) = x( uv ) ex P 



-M a u - -M b v 
2 2 



(3.34) 
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where x is an arbitrary function. The final result is 



exp 



M a r\ 



M b rl 



2 t n 



2 t n 



(3.35) 



We note that analogously to the conformal result pSfl , the scaling dimensions have 
to agree, while in this case we do not have a constraint on the masses A4 a ,b, since the 
system is not Galilei- invariant. 

The function x is partially determined from a few consistency conditions. We 
should expect to recover the bulk behaviour for large distances to the surfaces, that 
is for r a rb/r — > oo. Therefore, up to normalization 



X(u) ~ S MaMb e 



M a U 



u — > oo 



(3.36) 



On the other hand, for a free surface where the field vanishes, we expect absence of 
any correlations and thus 

X (0) = (3.37) 

Indeed, this is exactly the behaviour obtained from the method of images. We 
have for the surface correlation G s in terms of bulk correlations Gb 



G s (r,r';r) = G b {r - r',r) - G b {r + r',r) 
Gqt exp 



exp 



M (r + r 



A2 



_ / rr \ 
2G t x sinh I M — I exp 



Air 2 + r 



1 2 



2 T (3 - 38) 
and we identify x( u ) = 2G sinh(A^u) in agreement with the consistency conditions 
eqs. (p.36U3.37|) . 

Finally, if we were to impose in addition translation invariance in space, invariance 
under Galilei transformations is also implied and we do recover the bulk result ( [3.12|) 
for the two-point function. 



3.4 Two-point function for a non-stationary state 

We now consider a situation with a boundary condition at a fixed time. Boundary 
conditions of this type will be kept invariant by the subalgebra Sy together with the 
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scale transformation X , or rather its finite-dimensional subalgebra {X ,Y_ 1 / 2 ,Y 1 / 2 } 
for quasiprimary fields. For example, this may correspond to the situation when a 
system is in a predefined initial state and relaxes towards its critical equilibrium state 
0. Consider the two-point function of quasiprimary fields 

F = F(r a , r b - t a , t b ) = (cf> a (r a , t„)$(r 6 , t b )) (3.39) 

Invariance under space translations implies F = F(r; t a , t b ) with r = r a — r b . We next 

demand invariance under Galilei transformations, with e a = —e b = 1 

b 

J2 {UDi + eMin) F 

i=a 

= {{t a - t b )d r + M a r a - M b r b ) F = (3.40) 

In analogy to what was done before, this implies the Bargmann superselection rule 
M a = M b and 

F(r;t a ,t b ) =G{t a ,t b )ex P (-^— -J (3.41) 
Scale invariance demands that 

E (tidi + ^r,A + \x?) F = (3.42) 



Inserting ( |3.41| ), we find 

G{t a ,t b )=t-^ + ^/ 2 <$>{t a /t b ) (3.43) 

where $(t>) is an undetermined function. The final result is 

F = 5 MaMb r a ^^{t a /t b ) exp [-^ (J J^) (3-44) 

Note that here we have no condition on the exponents because the system is not 
invariant under the special Schrodinger transformation X\. 

A few consistency conditions may be added. If both times t a , t b are very large, we 
should expect to recover the form eq. ( |3.12| ) for the fully infinite space-time. In that 
case t a /t b — > 1 and we should have 

Hv)~6 Xa , Xb (l-v-iy iXa+Xb)/2 , v^l (3.45) 
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On the other hand, if tf, is kept finite, one should remain in the intermediate time 
scalime regime, termed "initial critical slip" in 1. In this case, one should expect 



~ v @ , v -»• oo (3.46) 

where 6 is related to the independent "slip exponent" defined in || . We shall present 
an example of this in section 4. 

3.5 Ward identity 

We now consider the effect of arbitrary coordinate transformations on correlations. 
We suppose that the system under consideration is described by a local action in d+1 
dimensions. This is obviously satisfied for static, but strongly anisotropic systems 
with local interactions. For many dynamical problems in d space dimensions which 
are at first defined via their equation of motion, there exists an equivalent equilibrium 
problem in d + 1 dimensions, usually supplemented with a "disorder conditions" to 
maintain the strong anisotropy fl31"| , |32[| . Then, considering the change of the action 
induced by an arbitrary coordinate transformation, the following identity holds, see 
e.g. § 

n 

^(0i(fi,ti) . . .64> a (f a ,t a ) . . .4> n (r n ,t n )) (3.47) 



a=l 



/ 



where we implicitly assume that in the correlators the Bargmann superselection rule 
is satisfied, the time T is denoted as the zeroth component of the coordinate R and 
Tij is the stress-energy tensor. As eq. (|3.47|) is written, we assume 8<p to contain 



all variations of the field (including changes of its phase) and to describe all 
changes of the action used to calculate the averages (...). The discussion presented 
here remains at the formal level of the equations of motion. We discard the possibility 
of anomalies which may arise from renormalization effects. Detailed discussions of 
these are available for conformal invariant theories, see Q, @], but for Schrodinger 
invariance, the analogous developments have not yet been done. 
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Since correlations are supposed to be invariant under infinitesimal Schrodinger 
transformations, we obtain a few constraints on the stress-energy tensor in complete 
analogy with the corresponding results for conformal invariance 01. The form of the 



8fj is taken from the generators eq. ( |2.6|) . Rotation invariance implies that is 
symmetric in space 

Tij = Tji , i,j = 1, . . . , d (3.48) 
Scale invariance gives the "trace condition" (here written in Euclidean form) 

d 

2T 00 + £T ii = (3.49) 
i=i 

which is the analogue of the vanishing trace condition in conformal invariance and 
the factor 2 comes from 6 = 2. Eq. ( [3.49 ) is satisfied in free non-relativistic field 



theory |Tj| (where it is written in Minkowskian form). Interacting non-relativistic 
field theories may give rise to anomalies, see |B3J. From Galilei invariance we find 



T 0i = , i = l,...,d (3.50) 

The requirement of special Schrodinger invariance does not add any further condition 

on T^. We have thus seen that 

translation invariance in space and time 
rotation invariance in space 
anisotropic scale invariance with 6 = 2 
Galilei invariance 
local interactions 



Schrodinger invariance (3.51) 



in analogy to the conformal result ||. In fact, using formally the equations of motion, 
we may even verify invariance under the entire infinite algebra S. This is for the time 
being the only indication that the generalization beyond iSfi n might be sensible. Again, 
the same type of result also holds for conformal invariance (when the central charge 
vanishes) . 

We do not go into a discussion of the possible anomaly structure here. As a 
preliminary exercise to that, we show in Appendix B that the Schrodinger algebra 
eq. (|2.7|) does not admit any non-conventional central extension besides the familiar 
Virasoro form for the generators X n . 
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3.6 Summary 

The main results of this section are the explicit expressions for the two- or three-point 
Schrddinger-covariant correlations in either an infinite or a semi-infinite geometry as 
given in eqs. (|3 . 1 . 35|,|3~4lD . Although the general form is quite similar to the 



corresponding results found from conformal invariance |13|, |28], |29[, there are some 
properties which come from the non-relativistic nature of the symmetry. The first 
one is the Bargmann superselection rule [EBJ for the masses. Secondly, the space 
dimension d has mainly the role of a parameter, at least for the quasiprimary fields 
only considered here, whereas the non-trivial group structure capable of extension to 
an infinite-dimensional algebra only occurs in the "time" coordinate. 

We remark that there is a certain analogy between Schrodinger invariance and 
conformal invariance close to a free surface j2l| . In both cases, the pair of com- 
plex linear projective transformations characteristic for full conformal invariance gets 
replaced by the subgroup of a single real linear projective transformation. 

The results obtained only use the finite-dimensional algebra. It remains an open 
problem how to extend Schrodinger invariance to the full infinite-dimensional algebra 
and find the scaling functions. 

4 Tests of Schrodinger invariance in exactly solv- 
able models 

We now turn to test the predictions obtained in the last section for some correlations 
in the context of some exactly solvable strongly anisotropic critical systems. We do 
not include here among the tests the well-known fact that the Green's functions of 
the free Schrodinger equation and the diffusion equation reproduce eq. ( |3.12| ) in any 
space dimension with x = d/2. We also refrain from discussing the range of possible 
applications of the models considered. 
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4.1 Kinetic Ising model with Glauber dynamics 

Consider the time-dependent ID Ising model with the classical spin Hamiltonian 7i = 
—JJ2i=i SiSi + i and Sj = ±1. To describe the time dependence, following Glauber |j34 | 



consider the probability distribution function P{s\, . . . , sl', t). It is often convenient 
to describe the evolution of P in terms of a master equation (see [^7], ^ for a detailed 
discussion) 

^P( 8l ,...,8 L ;t) = (4.1) 



(^2 w i( s i)^j P ( s l> • • • ) s L'i t) + Wi(-Si)P(s u 



• • • j Si, i) 



where the W{ are the rates describing the transitions between spin configurations. The 
following consistency conditions have to be kept. The first is probability conservation 
when summing over all configurations {s} = (si,...,s.l) 

T, p ({s};t) = i (4.2) 

{s} 

to be kept at all times. Second, the equilibrium distribution P eq ~ e _/3W has to be a 
stationary solution of the master equation, where (3 is the inverse temperature. This 
is usually implemented via detailed balance 

Wi(-Si) exp[-/3Jsi(si_i + s i+1 )) 



Wi(si) exp[/3Jsj(si_i + s i+ i)] 
Finally, averages are obtained from 



(4.3) 



<X>(t) = J2X({s})P({s};t) (4.4) 

{4 



Glauber |34| showed that the particular choice 



ro , (s , Hf(l _i^ 1+s , +l)) ( , 5) 

where a is the constant transition rate and 7 = tanh(2/5J), renders the model com- 
pletely integrable. We are interested here in his result for the time-delayed (con- 
nected) two-point function when the system is in thermal equilibrium at temperature 
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T \M 



G{r 1 -r 2 ,t 1 -t 2 ) =< s Pl (*iK(t 2 ) > c = e- at J2v lr - il Ie(ait) (4.6) 

t 

where r = r\ — r 2 , t = ti — t 2 , r\ = tanh(/?J), 1% is a modified Bessel function and the 
sum extends over the whole lattice. If there are no correlations between spins in the 
initial state, only the term with £ = r would be present. 



To analyse this, recall the asymptotic form, as x — > oo [36 



I t {x) ~ (2ttx)- 1/2 exp [x - (l + 0(x~ x )) (4.7) 
and obtain in the scaling limit r — ► oo, t — > oo with u = r 2 /t fixed 

G(r,t) * e-^)»'(2, lat )-^{expf-^) + g^exp(-^±|! 




~ (27rat)" 1/2 exp — — ) (4.i 

where in the last equation we performed the zero-temperature limit (since the ID 
Ising model has its critical point at T = 0). This holds exactly for vanishing corre- 
lations in the initial state and up to scaling corrections otherwise. This is indeed in 
agreement with the predicted two-point function eq. ( |3.12 ) and we identify x — 1/2 
and M = l/(2a). 

4.2 Lattice diffusion with exclusion 

Consider a system of many particles performing random walks on a chain of L sites. 



Each site can be either empty or occupied. The dynamics is defined as follows p2 |. 
First, pair all neighboring sites which can be done in two ways to be labelled A, 
chosen at odd times and B which is chosen at even times. At every time step, the 
dynamics of each pair is as follows. If both sites of a pair are either occupied or empty, 
the state of the pair is unchanged. If one site is occupied and one empty, the particle 
moves to the empty site with probability p or stays where it is with probability 1 — p. 
This stochastic rule for updating is applied in parallel to all pairs. 
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The time-delayed particle-particle correlation is 



G(r,t) =< n(r,t)n(0,0) > c 



(4.9) 



where n(r, t) = 1 if the site r is occupied at time t and n(r, t) = otherwise. The 
average is first over realizations of all possible time developments and second over 
an ensemble of initial states that is stationary with respect to the stochastic process 
32]. G was calculated exactly first for p = 1/2 |32]| and later for p arbitrary |37j with 



the result, for even times t and even sites r (similar results are known for the other 
cases, see p7[) 

(t-M)/2 



G(r,t)=p(l-p) 



+ E 



k=l 



(t-r)/2 
k 



(*-2 + r)/2 ) p t-2 k(1 



k-1 



p) 



21: 



(4.10) 



where p = N/ L is the particle density and N is the number of particles on the lattice. 
The limit L — > oo such that p is kept fixed is understood. In the scaling limit r, t — ► oo 



such that m = r 2 /t is kept fixed, this simplifies to |32|, ^7 | 



G(r,t) = p(l -p)(2 7 r J Dt)- 1/2 exp 



2Dt 



(4.11) 



and the diffusion constant D = pj (1—p). This agrees with the Schrodinger invariance 
expectation eq. (|3.12|) , and we have x — 1/2 and M. = 1/D. Similar results were 
obtained when quenched weak disorder is added, and it was shown that if the spatial 
disorder correlations decay rapidly enough, the same scaling form (|4.11|) results with 



a modified value of the diffusion constant D, see [37] for details. 

These results werde derived [J7j by mapping the system onto a six-vertex 



model satisfying a disorder condition. From the transfer matrix T = the quan- 
tum Hamiltonian can be obtained. In the time continuum limit r — > oo such that 
pr remains constant, the quantum Hamiltonian is found []37j to agree with the quan- 
tum Hamiltonian obtained directly from the master equation written in the form 

d t P = —HP, and reads 
/. 

2 



H = -\il [^X +1 + afaf +1 + AaX+i + (1 - A) (*? + <h) + A - 2] (4.12) 



24 



where the a x,y,z are Pauli matrices and with A = 1 for lattice diffusion. We note 



that this Hamiltonian, but now with A = 0, is also obtained [|3£| |39| from the master 
equation of the kinetic Ising model at T = considered above. The disorder condition 
for the vertex model formulation of lattice diffusion above makes the model undergo 



a transition of Pokrovsky-Talapov type [4~0|. In view of the common scaling form for 
both A = or 1, and because the spectrum of H is known to be A— independent 
we should expect this scaling form to hold independently of the value of A. 



4.3 Lifshitz point in the spherical model 

We now consider a static, but strongly anisotropic system. The model is the an- 
isotropic next-to nearest neighbor spherical (ANNNS) model, see [[|T| and references 
therein. Conventionally, it is defined by the Hamiltonian, on a hypercubic lattice 

Hsm = -E^-^+/ 5 " 1 cE^ 2 ( 4 - 13 ) 

i,j i 

where a, are real numbers, the spherical parameter ( is determined from the constraint 
< J2i s 1 >= A/" where M is the number of sites. Consider the model in D — d' + d 
dimensions. The couplings J?f are defined as follows. First, in all D dimensions, 
there is a ferromagnetic nearest neighbor interaction of energy J > 0. Second, in the 
d "parallel" directions, there is along the axes an interaction of energy kJ between 
next-to-nearest neighbors. The Fourier transform of J^f is 

J(k) = 2 J cos ^ + k E cos 2^") (4.14) 

\j=l i=l ) 

If a phase transition occurs, there is at k = k c = — 1/4 a meeting point between 
a paramagnetic, a ferromagnetic and an ordered incommensurate phase, which is 
referred to as a Lifshitz point 0, which is a strongly anisotropic critical point. When 
d' = 1, one has 8 = 2 and those coordinates in the d "parallel" directions will be 
referred to as "space" , while the one remaining direction will be referred to as "time" . 
In order to make the Galilei invariance explicit, we consider here a variant of this 
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model, which gives the same thermodynamics. The Hamiltonian is 



(4.15) 



where now s» is complex and the spherical constraint is 

< \s?\ 2 >= 2AT 

i 

We introduce the Fourier transform of the spin variable 



(4.16) 



s s = (2n)- D / 2 / dk H j l - n 



and get 



U 



dk 



- / dkA(k) 



The partition function is Z = J T>se l3n . Since 



< n % >=< n\ >= , <|/^| 2 >=A \k) 



(4.17) 



(4.18) 



(4.19) 



the spherical constraint ( ]4.16| ) becomes / dkAr l (k) = 2J\f, which is exactly the same 
as obtained from TCsm ( see e -g- EH) an d the free energy is F = — \nZ = 2F SM , 
where F$m is the free energy obtained from the Hamiltonian TisM- Consequently, 
the model defined by Ti is in the same universality class as the one defined by H-sm- 
In particular, the critical point is characterized by the condition ( = f3JD. The lower 
critical dimension is at D = 3, the upper critical dimension at D — 7, provided d' — 1 
which we assume from now on. 

Consider the two-point function C(a — b) = 9ft < s^s*- > where 9ft denotes the real 
part and 

k B T 



Z(2tt) d J ~ *~ J ^k^t 
(2n)- D k B T I dkh7 l {k)j t{S - %) 



(4.20) 

This has been calculated for integer dimensions D [[i3[| . In the scaling limit r — > 
oo,t — > oo with r 2 /t fixed (r and i are the distances in space and time between the 
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points a and b) the result is 

C{a-b) = C(r,t) 

A 2 = 

The function \P(a, x) is given in table p] 



^2 r (D-3)/2^ 



J 



2D-3 



7T 



JJ-3 r 2 ^ 
4 '27 J 

1 \ n D-2 

4 



(4.21) 



a 


ty(a,x) 


asymptotics 


l 

4 
1 
2 
3 
4 

1 


1/T(3/ 


I-i/i(x/2) + J 1/4 
t 2 x/2) 1 / 4 [/_ 1/4 (x 

'4) + v^F(x/2) 3 / 4 


(x/2) 
)-L 

L l/4( 


_i/ 4 (x) 

ar) - ii/4(z) 




2X" 1 / 2 
(2/r(l/4))x- 1 

-l/(2r(3/4))x" 2 



Table 1: Scaling function \l/(a, x) at the Lifshitz point of the spherical model and their 
leading asymptotic behaviour for x — > oo. I n and K n are modified Bessel functions 
and L n is a modified Struve function. 



Comparison with Schrodinger invariance eq. ( |3.12| ) shows agreement for the case 
D = 6, while the scaling function has a different form in the other cases. Recall that 
the system given by the Hamiltonian Ti has a dispersion relation of the form 



E — 



Am? 



E + —k 2 
2m 



E-—k 2 
2m 







(4.22) 



rather than E = k 2 /(2m) which was used for making the Schrodinger-invariance 
predictions. We can only expect to recover the results of Schrodinger invariance if 
the propagator actually solves the free Schrodinger equation, and not just a more 
general forth-order differential equation. We see this to be the case for D = 6 from 
table [I] and concentrate on this from now on. 

In order to test the prediction for two- and three-point correlations, we consider 
several scaling fields as defined in table 0. We also give the scaling dimensions and 
the masses (in units of the mass of the field a) of these fields. Concerning the field 
a, we can confirm its values for x a and v a from eq. ( 4.21J) . 
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Oa 


Sa 


3/2 


1 


A 


£a 




3 


2 


V2A 2 


Va 




3 





A 2 


Sa 




9/2 


3 


Vqa 3 



Table 2: Some scalar scaling fields arising at the Lifshitz point of the spherical model 
at D = 6 and their scaling dimensions x^, masses M.^ in units of the mass hA a and 
the normalization such that v~^ x <\> has $ = 1 in eq. ( [3.121) . a' denotes a space-time 
point on a neighboring site of the point given by a. 

— * — * 

The calculation of the other correlations is simple because, since A(— k) = A(/c), 
the imaginary part of < s a s* b > vanishes and we have C{a — b) =< s^si >. 

Consider the two-point functions first. (From now on, a, b, c denote space-time 
vectors and we also write a = (r a ,t a ) etc.) Obviously, < e a >= 0. Then, 



<e a e* b > = < s a s a ,s* b s* b , >= C(a-b)C(a' - &') +C(a-b')C(a' -b) 



2[C(a-b)] 2 = 2A 4 t- 2x °e- r2/t 



(4.23) 



where in the second line the scaling limit was taken and we confirm the result given 
in table 0. Next, we have < r] a >= C{a — a') = const. Then the connected two-point 
function is 

< VaVb >c=< s* a s a ,s* b s bl > c = C(a - b')C(b - a') ~ t~ 2x °A 4 (4.24) 
whereas the exponential terms cancel. Finally, we have < S a >= and 

< £ a ££ >=< s a s a ,s ar/ s* b s* b rsl„ >~ 6[C(a - b)] 3 = 6A 6 r 3x ° e - (3/2)r ' !/t (4.25) 
and we have verified all entries in table [| It is straightforward to verify that 

< a a e* b >=< a a 7] b >=< aj^l >=< e a T] b >=< e a ££ >=< 7] a Z* b >= (4.26) 



This fully confirms the Schrodinger two-point function eq. ( |3.12| ). 

At this stage, we clearly recognize the difference between the models described 
by Ti-sM and 7i, respectively, from the point of view of Schrodinger (in fact already 



28 



Galilean) invariance. The two fields e and rj have the same scaling dimension, but 
different masses. In the context of the model described by TCsm, however, these 
two distinct fields get lumped together into just one field, e, say. Correlations of e 
will in general not satisfy Galilean invariance. On the other hand, there is just a 
single order parameter field a. Since in many applications one is only interested in 
the order parameter correlations with itself, it is for that restricted purpose enough 
to stay with the conventional form of TCsm, rather than go to 7i with the correct 
Galilean transformation properties. 

We now turn to the three-point functions. To check the Bargmann superselection 
rule, one may verify that, for example 

< s a s h s* c >=< e a e b e* >=< i] a a b e* c > c =< cr a cr b £* >=< e a e fe S* >= (4.27) 

A non- vanishing correlation is 

< cx a a b e* c >=< s a s b s* c s* c , >= C(a-c)C(b-c r )+C(a-c')C(b-c) ~ 2C{a-c)C{b-c) 

(4.28) 

Since x e = 2x a , this agrees indeed with the prediction eq. ( |3.28| ) from Schrodinger 



invariance and we identify the scaling function ^ aa ^ = a/2 which is a constant and 
we have used the normalization given in table |2|. Another check of the Bargmann 
superselection rule is provided by showing that < e a a b cr* >= 0. Furthermore 

< a a e b E* c > = < s a s b s b ,s* c s* c ,s* c „ >~ QC(a - c)[C(b - c)f (4.29) 

= QA 6 (t a — t c )~ Xa {t b — t c )~ 2X,J e~ ( - Ta ~ Tc ^ 7 2 (*a-*e) e -(n-rcf '/(fb-tc) 



which from table ^| is seen to agree with ( 3.28|) . For normalized fields we identify 



^cre,E = \/3- Apparently, the massive scaling fields do reproduce the predictions of 
Schrodinger invariance. 

Finally, we look at some examples with the massless field i] whose correlations are 
not immediately zero. For example 

< VaO-to-* >c=< s* a s a >s b s* c > c ~ C(a - b)C(c - a) (4.30) 
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We denote a = (r„, t a ) and get 



< VaW* >c= A\t a - t b )- x "{t a - t c )-^ e - (ra - r6)2/2( * a - t6) e +(ra - rc)2/2( *°-* c) (4.31) 



which indeed agrees with ( |3.28| ) (recall from eq. ( j2.8| ) that a* picks up a phase opposite 
to cr), since x v = 2x a . Using normalized fields, we identify ^ V a,a = 1. Next, consider 



< r] a e b e c > c 



< s* a s a /SbSb>s* c s* c , > c ~ AC(a — b)C(c — a)C(c — b) 
AA 6 [(t a - t b )(t a - t c )(t b - t c )]- x ° exp 

.„-(r c -r b ) 2 /2(t c -t b ) 



(4.32) 



(r a - r b ) 2 (r a - r c ) 2 



2(£a ~ tb) 2(t a — t c 



To relate this to the form ( p.28| ) , we note that the first line of the result fl4.32| ) already 
takes the expected form and we merely have to rewrite the last factor in the required 
scaling form. Consider the argument of the scaling function \1/ 

[(r c - r a ) (t b - t a ) - (r b - r ) (t c - t a )} 2 (r c - r b f (r c - r a ) 2 (r b - r a ) 2 



(tb ~ t a )(tc — t a )(tc ~ tb) 



tr — th 



■ + 



tr t n 



th — t„ 



The last two terms on the right-hand side can be rewritten as 



(r e - r a ) 2 (r b - r a ) 2 (r c - r a f 



tc t a t b t a t c t a \ \?"c TfiJ \t b t a 



'"ft ^ a \ ( tc t a 



(4.33) 



(4.34) 



in the scaling limit. We thus find agreement with the expected form and identify for 



the normalized fields ^f vei£ (u) 



2exp(— u/2). Finally, we consider 



< VaVbVc >c = < S a S a 'S b S b 'S c S c > > c 

~ C(a-b)C(c-a)C(b-c)+C{a-c)C{b-a)C{c-b) (4.35) 

= A 6 [(t a - t b ) (t c - t a ) (t b - t c )r X ° ■ 

{exp f-~[A , 6 + A c , a + A 6 , c ]) + exp f-~[A 6>a + A a , c + A c>6 ]) J 

where A a6 = (r a — r b ) 2 / (t a — t b ) = — A fea . The very fact that this correlation does 
not vanish again confirms that the field r] is massless. We verify that 

A 6 + A + A b K r " ~ r c)(tb - t c ) - (n - r c )(t a - t c )] 2 ^ 3g ^ 

(t a — t b )(t a — t c )(t b — t c ) 

and find, identifying ty mjT) (u) = 2cosh(w/2), complete agreement with the prediction 

eq. (p|). 
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4.4 Critical relaxation of the spherical model 

Having studied in some detail the case of a setting of infinite extent in both time 
and space directions, we now look into one example with a macroscopically prepared 
initial state which is not the equilibrium state and the system is then allowed to 
relax towards equilibrium. While for very long times, we are back to the dynamical 
scaling considered so far, it was realized || that already the intermediate stages 
of the relaxation process display universal behaviour. The exponents and scaling 
functions describing this have been calculated via an e— expansion ||. Since we are 
merely interested in the special case of a dynamical exponent z = 2, we concentrate 
on the n — > oo limit of the 0(n) vector model. It was shown in |J that for a 
purely relaxational dynamics towards the equilibrium critical point the exact response 
function is in d < 4 space dimensions 

G^t, t') = Q(t - t') (t/t') 1 ~ d/4 exp [~Xk 2 (t - t')] (4.37) 

where Q(t — if) is 1 for t > if and zero otherwise. Fourier transformation in space 
gives 



G(r ] t,t>) = Q(t-t')(n/\) d / 2 r d / 2 (1) 



t \ ^ ( t , 

1- T exp 



d/2 

(4.38) 



4A 2 t - V 

This is indeed consistent with the result eq. ( |3.44 ) for the two-point function. We 
identify x = d/2, M = 1/(2A 2 ) and 

$(v) = e(v - l)(ir\) d ^v 1 - d ^ (l - v- i y d ' 2 (4.39) 

which also agrees with the consistency conditions given in section 3. 

5 Some remarks beyond 6 = 2 

Since for generic values of the anisotropy exponent 9 there is at present no general 
approach available, we content ourselves with a few results from selected models. In 



any case, these permit to submit the conformal invariance result eq. ( |1.4|) [ 10[ to a 
test. 
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5.1 Lifshitz points of higher order in the spherical model 

We come back to the ANNNS model introduced earlier. Now, we add further in- 
teraction terms along the axes denoted as "space" dimensions. Since we are only 
interested here in the spin-spin correlation, it is sufficient for us to consider the real 
Hamiltonian TisM eq. Q4.13Q - The Fourier transform of the couplings now is 

(D d n \ 

5^cosJfei + KiCos[(i + l)kj}\ (5.1) 

i=l j=l i=l J 

Previously, we had taken k = k% to be the only non-vanishing coupling. With several 
of the Ki non-zero, the phase diagram will contain lines of Lifshitz points (also called 
Lifshitz points of first order, where L = 2) which end in a Lifshitz point of second 



order (with L = 3), in analogy to the definition of multicritical points, see [41] for a 
review. Lifshitz points of higher order are defined analogously. At a Lifshitz point of 
order L — 1, we have 

n i D d 

J(£)~2J£> + E fci"CiX>f + ••• ( 5 - 2 ) 

i=l 1 j=d+l j=l 

which defines the readily calculable constant cl- When d! — 1, the anisotropy ex- 
ponent is 9 — L. We are interested in the critical correlation function C(a — b) = 
< a a (Jb >■ This can be calculated exactly, see f43|. As we had already seen for 



Lifshitz points of first order above, the correlations of the model considered here will 
only in special cases solve the dispersion relation E ~ k e , rather than E 2 ~ k 2e . This 
should be the case if the scaling function does not just show a power-law behaviour 
for large values of u. This is the case ii|] 

D = L + 2 + 2m,m = l,...,L-l (5.3) 

Then the critical two-point function is in the scaling limit r — > oo,t — > oo with 
u = r e /t fixed |3| 



Cif A = A L ,^l^[L^^ L ^ (5.4) 



3 Eq. (2.22) in contains a typing error and correctly reads a=\d— -377^ — 1 = |(d — d_) 
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where H can be expressed as a finite sum of generalized hypergeometric functions and 
Al,d is a known non-universal constant. Here we merely consider the behaviour for 
large values of the scaling variable u, to leading order 



C(f,t) ~ BL ^ t -frn+l)/L u {D-3S)/M6-l)) 

exp 



71 L \ „i/(9-i) 



cos 



CiCos( 2L3ir 

sin I «V(9-D a. ![ L D ~ 3L 



(5.5) 



.2L-17 2L-1 2L 

where and > are known non-universal constants. We remark that only the 
power prefactor of the scaling function depends on D, the rest of the scaling function 
depends in the large u limit only on L. 

Let us compare this result with the prediction eq. ( |1.4|) following from confor- 
mal invariance flQfl . While conformal invariance gives for a non-conserved order 
parameter a simple exponential behaviour for the scaling function, we rather find 
for the higher order Lifshitz points (L > 3) a stretched-exponential behaviour. The 
two forms only agree for 9 = 2. On the other hand, for a conserved order pa- 



rameter the conformal invariance scaling function for u large is of the form |T0 
ip{u) ~ ^-sz/s^g-u 1 / 3 cos^y^^i/S) w hj c h for z = 4 is van Hove theory. It is interesting 
to note that for L = 4, the ANNNS model reproduces the same behaviour. 

Finally, we look into the case where 9 = 1/2. This is realized if d = 1 and 
L = 2. Now the direction parallel to the next-to- nearest neighbors interaction will be 
interpreted as "time" and the other direction are referred to as "space". For D = 4, 
we find again an exponential-like behaviour for the spin-spin two-point function |^3 

C(f, t) ~ r^ 2 u- 3 exp (-I«V(*-U) (5.6) 
in the scaling limit with u = r 1//2 ft fixed and 9 = 1/2. Again, this is different from 



the conformal invariance prediction eq. (|1.4|). 
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5.2 Directed percolation 

As a further example for a strongly anisotropic critical system we consider directed 
percolation in 1 + 1 dimensions (for a review, see ||). In percolation, sites or bonds 
are filled at random with probability p and percolation proceeds along paths between 
occupied nearest neighbor links. In directed percolation, there is in addition a pre- 
ferred direction and percolation is allowed to proceed only in one sense along this 
direction. This preferred direction is called "time" and the orthogonal ones "space" . 
Consider the pair connectedness G(f,r') which is a measure of the probability that 
sites at r and r ' are connected by a percolating path. It is well known that there is 
a critical value p c such that one has the scaling form 

G(r,t) = Ar 2f} l v W®{Bv) , v = r/t^ (5.7) 

where (3,v±,v\\ are critical exponents, the anisotropy exponent 9 = v\\jv±_ and r,t 
measure the "space" and "time" distances, respectively. Precise numerical values on 



p c on various lattice and of the exponent have been obtained, see H3, M. For our 



purposes it is enough to notice that in (1 + 1)D, to which we restrict ourselves here, 



we have 9 ~ 1.5807 |5 . 



Benzoni [4"1] studied the pair connectedness by calculating numerically the mo- 
ments 

x (») = A r dv\v\ n $(Bv) (5.8) 

J — oo 

and verified that the following ratios 

[ Y (1)12 [ y (2)l2 r (2)12 

y(0)y(2) ' r y(0) y (4) ' ^ y (l)y(3) 

A A A A A A 

are independent of the non-universal scale factors A, B and should therefore be uni- 
versal. A careful numerical computation 0] then yields numerical values for C, F, G 



for various lattices of both directed site and directed bond percolation. The results 
are in full agreement with universality p4 |. 



We proceed to analyse these results in the following way. We try the ansatz for 
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the scaling function $(t>) 



$(v) 



\v\ exp(— v 



(5.10) 



where a, b are constants to be determined. Then x — 2a + b + We 



now fit this form to Benzoni's MM numerical results for C, F and G and find 



a = 2.49 ±0.16, 6 = -0.016 ± 0.03 from C 
a = 2.58 ±0.12, 6 = -0.023 ± 0.03 from F 
a = 2.62 ± 0.17 , 6 = -0.023 ± 0.07 from G 



(5.11) 



with the mean values a = 2.56(7) and b = —0.02(3). Since the scaling function $(t>) 
is finite for v = we interpret this result as implying that b = 0. In table |3| 
we give the results for a as found using the ansatz (|5.10|) with 6 = from the ratios 



ratio 


1 


2 


3 


4 


5 


C 


2.41 


2.56 


2.49 


2.24 


2.39 


F 


2.50 


2.69 


2.69 


2.32 


2.49 


G 


2.54 


2.74 


2.65 


2.37 


2.54 



Table 3: Exponent a as determined from the moment ratios C, F, G for the following 
realizations of directed percolation: (1) square bond, (2) square site, (3) square site- 
bond, (4) triangular bond, (5) triangular site. 

C, F, G and various realizations of directed percolation. Note that the estimates 
for a obtained from different moments and different lattice realizations of directed 
percolation are the same, which means that the chosen ansatz does indeed describe 
the available data. From all this we conclude 



a = 2.6 ±0.2 



(5.12) 



Making contact with our previous results, formulated in terms of the scaling variable 
u = v 6 = r 9 /t, we obtain from ( 5.10|) 





a 



9-1 



~ 2.72. . . 



(5.13) 
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using the known value of 9. Comparison of ( 5.12Q and ( 5.13Q implies that also in this 



class of models the two-point function scaling function appears to be consistent with 
the same stretched-exponential form as already observed for the Lifshitz points of 
the spherical model and in disagreement with conformal invariance eq. ( |1.4| ). 

6 Conclusions 

In this paper, we have examined the simplest consequences of the hypothesis of 
local dynamical scaling with space-time dependent local rescaling factors A(r, t). We 
have seen that for the special case of an anisotropy (or dynamical) exponent 6 = 2, 
the Schrodinger group, which is the non-relativistic limit of the conformal group, is 
a sensible candidate for a group of local scale transformations. The treatment of 
Schrodinger invariance (of quasiprimary fields) is in many respects quite analogous 
to conformal invariance. However, there are a few distinctions, the main one being 
the role of the phase transformation which is not present in the conformal group. We 
hope that the experience obtained in this simplest non-conformal case may become 
useful for the extension of the method to generic anisotropy exponents 6. We have 
derived the form of the two-point and three-point functions for both infinite space 
and time (eqs. (|3.12| , |3~28|) ) and for the two-point function also if either space or time 



is restricted to the half-infinite space (eqs. (|3.35|]3.44| )). The results obtained are in 
agreement with and extend those following from the weaker restrictions of Galilean 
invariance. 

The Lie algebra of the Schrodinger group can be naturally extended to an infinite- 
dimensional one. We have not solved the problem of how to use this infinite algebra 
to calculate the critical exponents and the scaling functions in the correlations which 
are left undermined in this work. We hope to come back to this in the future. 

Several exactly solvable statistical models with anisotropy exponent 6 = 2 were 
seen to reproduce the results of Schrodinger invariance for the two- and three-point 
functions. In particular, we have seen that due attention must be paid for correctly 



36 



implementing the changes of the phases of scaling fields as demanded by Galilean 
invariance. 

Evidence from some models with anisotropy exponent 9^2 suggests that, at 
least for large values of the scaling variable u = r 9 /t, the two-point scaling correlation 
function might behave as 

$(u) ~ exp (6.1) 

(where oscillating and power-law prefactors as well as non-universal scale factors 
were suppressed). We have found examples for 9 = n, with any integer n > 3, for 
9 ~ 1.58 . . . and for 9 = 1/2 (for which the exact scaling function eq. ( |5.6| ) is known). 
This finding is in disagreement with the form eq. (|1.4j) suggested by using conformal 
invariance in space. 
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Appendix A. Solution of a system of linear differ- 
ential equations 



We derive the general solution H = H(r, s;r,a) of the following system of differential 
equations 

(rd T + ad a + ^rd r + ^sd^H = 

(rd r + ad s )H = (A.l) 
T 2 d T + a 2 d a + rrd r + asd s )H = 



The technique consists of subsequent solution and resubstitution, see J3"0| . The second 
of the eqs. (|A.1| ) is solved by 

H = H(t;r,a) , t = - - - (A.2) 

r a 

while the other two equations become 

(rd T + ad a -hd^H = 

(r 2 d T + a 2 d a ) H = (A.3) 

The first of those is solved by 

H = H(u, v) , u = rt 2 , v = at 2 (A.4) 
and the second one becomes 

(u 2 d u + v 2 d v ) H = (A.5) 

with the solution H = ^(u^ 1 — v~ x ) where \& is an arbitrary function. Backsubstitu- 
tion then yields the result eq. ( |3.27| ) in the text. 
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Appendix B. Impossibility of non- conventional cen- 
tral extensions 



Consider the centrally extended (infinite) Schrodinger algebra 



[Xn, X„ 
_X n , Y n 

[X n , M v 
\Y n , Y n 

\Yfij M n 
[M n , M n 



[n 



n 



- m)X n+m + ^- (n 3 - nj 5 n+mfi 



- - m ) Y n+m + D(n, m) 



—mM n+m + E(n, m) 
(n - m)M n+m + A(n, m) 
F(n, m) 
Kn5 n+mfi 



(b.i; 



where A, D, E, F are numbers and c, K are constants. The special form of the central 
extensions for X n , M n is well known. We show that, with the only exception of c, these 
central extensions either have to vanish or can be reabsorbed into the generators. 

This follows from the Jacobi identities. Begin with D(n, k). Consider [X n , [X m , Y^]] 
and their cyclic permutations. This implies 



— - kj D(n, m + k) - f - - k) D(m, n + k) - (n - m)D{n + m,k) = (B.2) 



Let n = 0. It follows (besides £>(0, 0) = 0) 

1m — 2k 



Dim, k) = D(0, m + k) =: k) dim + k) 

2 m + k V 2 ' 



(B.3) 



which defines d(k). Turning to E(n,m), consider [X n , [X m , Me]] and their cyclic 
permutations. This implies 



— £E(n, m + £) + £E(m, n + €) — (n — m)E(n — m, . 







(B.4) 



Now take m = in ( |B.4| ) and then either £ = 0orn + £ = 0, implying E(n, 0) = 
for all n. We therefore write E(n,m) = mS(n,m). Inserting in (p.4j ) and taking 
m = 0, we find that £(n, 1) = £(0, n + £) + 5 n+ ^ e(n) + ^o^^)- With the definitions 
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8(0, n) = e(n) and 77(71) = —ne(n), we have E(n,m) = me(n + m) + i](n)5 n+mt o. 
Backsubstitution into ([B.4j) implies 



(n + m) (rj(n) — r](m)) = (n — m)r](n — m) (B.5) 

Let n = and get 77(71) = r/(— n). Let n + m = to see that i](2n) = for all n. 
Taking m = 2n, we get r)(n) = —~r](—n). Consequently, t](n) = for all n. Turning 
to A(m, k) = —A(k,m), consider [X n , [Y m ,Y k ]} and permutations. We get 

(m — k)E(n, m + k) + ^— — k^j A{n + k, m) — ^— — rri^j A{n + m, k) = (B.6) 

We use the result for E(n,m), the antisymmetry of A(m, fe), insert in ( p.6|) , put 
n = and divide by m + Ho get A(k, m) = (m — k)e(m + k) + a(/c)5fc +mi0 , with 
a(fc) = — a(— fc). Backsubstitution into ( |B.6| ) then implies (3n/2 + m)a(—m) = 
(n/2 — m)a(n + m). We now choose n = 2m and get a(m) = 0. To see that K = 0, 
consider [F„, [Y m , Mfe]] and its permutations to get (n — m)Kn5 n+m+ kfi = 0. Finally, 
we turn to F(m, k) and consider [X n , [Y m , Mk]] and permutations to obtain 

kF(m, n + k) - - rnj F(n + m, k) = (B.7) 

Put n = to find F(m,k) = f(m)5 m+ k t o- We now have to distinguish two cases, 
i) The index m of Y m is half-integer. Since the index k of My. is always integer, we 
directly have F(m, k) = 0. ii) The index m of Y m is integer. Backsubstitution into 
( |B.7|) then gives 

(n + m)f(m) + ^— — /(n + m) = (B.8) 

Let n + m = and find /(0) = 0. Then, let m = and get /(n) = -2/(0) = 0. 

Consequently, the only surviving terms are given by d(m) and e(n). These can be 
absorbed into the generators by defining M n = M n — e(n) and Y m = Y m — d(m). The 
only central term remaining is the one parametrized by c. This proves the assertion. 
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